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Abstract 

We investigate the efficiency of acceleration of charged dust particles by low-frequency 
Alfven waves in nonlinear approximation. We show that the longitudinal accelera- 
tion of dust particles is proportional to the square of the soliton amplitude 0(|6mP)i 
while the transversal acceleration is oi 0{\bm\)- In the conditions of the interstellar 
medium the resulting velocity of dust particles can reach 0.3 to 1 km s~^. 
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1 Introduction 

The presence of charged dust grains in space plasma may play crucial influence 
on the overall dynamics, as for instance, the breaking down of frozen-in-field 
conditions [1] and dissipation of magnetic flux [2] in the interstellar molecu- 
lar clouds. The dynamical importance of dust grains, and circulation of their 
mass (dust destruction and growth) in the interstellar medium is essentially 
determined by their kinetic temperature, and therefore understanding of how 
they are accelerated is of significant interest [31H] . In general, magnetohydro- 
dynamic (MHD) waves are known as an efficient heating source of plasma 
components. For instance, interaction of dispersive Alfven waves with plasma 
results in Joule heating of the electrons, generation of zonal flows and forma- 
tion of localized Alfven structures and Alfven vortices [5]. In a dusty plasma 
MHD waves, in particular, compressible fast MHD modes, are found to be ef- 
flcient in heating of translational degrees of freedom of dust particles through 
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gyro- resonant acceleration [3] . More recently ponderomotive forces of the shear 
Alfven waves are shown to be a powerfull mechanism of dust acceleration [6], 
such that in a short time a microsized dust particle of mass ~ 10~^^ g can 
reach a few km s~^ velocity - comparable to the phase velocity of the Alfven 
wave. In [6] a steady nonuniform field envelope with a fixed field gradient 
(caused, for instance, by a nonuniform plasma density distribution along the 
wave vector) is described. Contrary, in our analyses below the ponderomotive 
force is unsteady in the rest reference frame where the pulse envelope moves 
with the group velocity. In this approach we are able to derive the perturbation 
amplitudes in the reference frame connected with the pulse, for densities and 
longitudinal velocities of plasma components driven by the mutual action of 
Coulomb and ponderomotive forces. Therefore, the nonuniformity of the field 
amplitude and the corresponding plasma acceleration in the pulse follows in 
our approach from a self-consistent solution of MHD equations. 

In [6] they considered a relatively high-frequency limit, corresponding to fre- 
quencies much higher than the dust cyclotron frequency Ud- It might be ex- 
pected though that in the low-frequency end u <ti ujd the presence of dust 
can have an important influence on dynamics of Alfven waves. In particu- 
lar, one can think that if dust particles are accelerated, the energy gained 
by them from the Alfven wave can be sufficient to decrease its amplitude, 
which in turn can reduce the acceleration efficiency. In this paper we explore 
the efficiency of acceleration of dust grains by the nonlinear low-frequency 
Alfven waves, and show that the nonlinearity strongly restricts the ability of 
Alfven waves to transfer their energy to charged dust grains. Such solitons 
can be actually treated as elementary objects of MHD turbulence in the low- 
frequency range. We follow therefore recent discussion of interaction of such 
low-frequency Alfven solitons with the plasma [7118] . 

In the conditions of the interstellar plasma this low-frequency range, uj -C 
~ 3 X 10~^^ Hz, can be of great interest from the point of view of dust ac- 
celeration because a considerable energy fraction of the interstellar turbulence 
can be locked in low- frequency motions [9]. Of course, nonlinear interactions 
move low-frequency Alfven perturbations toward higher frequences. However, 
it seems natural to assume that in the interstellar environment a steady state 
spectrum of perturbations establishes on a relatively short time scale (defi- 
nitely shorter than 100 Myr - the crossing time for a spiral wave), so that 
perturbations are always present in the low-frequency end. 

We proceed as follows. In Section 2 we find the dispersion relation and the 
group velocities for longitudinal Alfven and fast MHD waves, in Section 3 
we derive the nonlinear wave equation; Section 4 contains estimates of the 
velocities of dust grains accelerated by the Alfven and fast MHD soliton waves; 
summary of the results is given in Section 5. 
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2 Group velocity 



In a fully ionized dusty plasma the dispersion relation for longitudinal Alfven 
and fast MHD modes can be written as 

iV' = ei ± 62, 



where N = kc/u, 
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a = e, z, d, LUpa and are the plasma and cyclotron frequencies of the ath 
component. In the low- frequency limit uj ujd equation ([T]) can be written in 
the form 
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where signs ± belong to the Alfven and fast MHD waves, respectively, VAi = 
Bo/^/4^Fpi is the Alfven speed. 

Assuming mdPdl ZdfriiPi ^ 1 and uj <^ pdOJdj P-, P = Pi + Pd, one can get from 
here 

= k\,l (l ± ^kVAd] , (4) 



where v\ = Bq/Atip, r^d = VA/^d- The group velocity and its derivative for 
Alfven waves are 

u = — = va \ l ^ TAdk , -j7 = —VArAd > 0, (5) 
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and for fast MHD mode 
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3 Nonlinear parabolic equation for low-frequency Alfven and fast 
MHD waves 



Following the procedure described in [TUfTT] one can derive the nonlinear 
Schrodinger equation for the envelope of the Alfven wave package in the form 

dEi u'd'^Ei u ^ ujlk^^ ( kv,k \ ^ 

+ 7717^2^^ — — hEi = 0, 7 



where 



U = ± ■ rADVA, 

Pi + Pd 



is the dispersion of the group velocity; sign "plus" belongs to the left-polarized, 
"minus" to the righ-polarized waves, Nk = ftk/nko is the normalized density 
perturbation of ki\i component. In the left-polarized wave ujk is positive for the 
electrons and dust particles, and negative for the ions, in the right-polarized 
Alfven mode is positive for the ions and negative for the electrons and dust 
particles. 

Assuming vti <^ t'A ^ one can neglect the electron inertia and write the 
equation of motion in the form 

dNe _ e d(j) ( d k d \ |Ep 



The equations of longitudinal motion of the ions and dust are 

dVz ^ e_d(^ ^ {—+ ^ ^\ |^|2 

dt M dz M'^uiui — uj) \dz uj uji — u dt j ' 

dwz Zed(f) ^ Z^e^ |^ ^ _l_ ^ ^'^ \ |^|2 



Last terms in the r.h.s. of equations fl9l)- f|TTl) describe the ponderomotive force 
acting on the particles by the moving wave. 

With using expansion over a small parameter uj/ui -C 1 one can transform 
equation of motion of the ions and dust f|TOl) - f|TTl) in the reference frame moving 
with the wave group velocity to the form 

dV. edct^ f + F 
— u—— = — , 12 
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After integration of ^ and ffT^ - f[T^ one can find tlie density perturbations 

as 
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Then the quasineutrahty condition 

UioNi = UeoN^ + ZudoNd, (14) 



gives the following solution for the Coulomb potential 
e0 Wi W2 



and the explicit solutions for the density perturbations 
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here c^riioTe/neoM is the square of the ion-sound dust speed. The longitudinal 
velocities are connected through the continuity equations with the density 
perturbations by the equations 

v^k = uNk. (17) 



With accounting of equations (fT6l) and (fT7|l in the last terms in equation ([7]) 
one can reduce it to the following form for the new variable h = B/Bq = 
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The Lighthill criterion for the modulational instabihty fulfils when the signs of 
u' and of the nonlinear term in (fTSl) coincide. Thus, for Alfven wave with n' > 
modulational unstability takes place whenever the conditions c^pa/p < < 
hold. In the interstellar plasma with pd/p ~ 0.01 both these conditions are 
valid. For the right-polarized Alfven waves with u' < modulational instabil- 
ity can develop only in one of the two cases: either > cj. or < c^pd/p, 
which hardly can hold in the conditions of the interstellar matter. Therefore, 
left-polarized Alfven waves with the amplitude higher than a critical value 
experience in the interstellar plasma modulational instability which results in 
formation of a sequence of Alfven solitons. In the limit u ^ Cg the free term 
in (1181) becomes singular, which means that the quasineutrality condition is 
not valid anymore. 

Using the standard procedure we represent the magnetic field as & = fo_Le**, 
where b± = b±{( — vt), $ = $(C — ^i^), what means that the envelope solution 
we are seeking for drifts slowly (i.e. v <^ u) in the reference frame moving 
with the group velocity, while the phase $ drifts with the velocity vi ^ u. 
Separating the real and imaginary parts of (ITSll we arrive at 



d'b^ v{v-2v,) 2a 
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where 
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The solution of ffTOD is 



b± = b^ch-' (^^) , (21) 



where A = ^Ju'/ab'^, and vi can be found as vi = v{l — au'b'^/v'^)/2. In 
general, the solution (!2T]) is determined by the two free parametes: b^ and 
V. It is readily seen that collisions of dust particles with ions and atoms are 
unimportant on the scales of interest. Indeed, the characteristic length of the 
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soliton (I2T!) for uj <^ Ud can be estimated as ^ rAdb^, or for typical 
parameters in the interstellar plasma ^ lO^^b^ cm, which on the lower 
end is much less than the drag free path of dust particles: £ ~ 3 x lO^^n"^ 
cm, where the mass ratio rrid/m ~ 10^*^ and mean grain radius a = O.l/xm are 
assumed. 



4 Longitudinal and transversal acceleration of dust particles 



In the laboratory reference frame longitudinal acceleration of dust particles 
is determined by the last equation in (fT6|) and (fT7|) : \Vdz\ — u\b\'^/2, so that 
within the validity of the approximation the amplitude of dust velocity on 
nonlinear Alfven wave can be of several percents of the Alfven speed, which 
in the conditions of the interstellar plasma can vary from ~ 3 km s~^ in 
molecular gas (p ~ 10~^^ g cm~^, B ~ 10/uG) to ~ 10 km s~^ in diffuse 
HI phasei (p ~ 10~^^ g cm~^, B ~ 3yuG). Therefore, reasonable conservative 
estimate of dust acceleration along magnetic lines in such environments may 
be of 0.3 to 1 km s~^, respectively, which however depends on exact value of 
the soliton amplitude. This is much lower than the estimate obtained in [B]. 
Similarly, the dust velocity in perpendicular direction due to acceleration by 
the soliton can be estimated as \Vd±\ ^ u;\b±\/k u\b\, which is therefore a 
factor of larger than acceleration in the longitudinal direction. 

The acceleration of dust grains by a regular low-frequency Alfven soliton 
{uj ^ ujd) is connected with the ponderomotive force as in [6], and differs 
qualitatively from the mechanism described in [3], where translational heat- 
ing of dust particles is associated with stochastic acceleration of dust particles 
by MHD waves with ~ cj^ through cyclotron resonant interactions. Indeed, 
it is readily seen that dust particle moves under the action of soliton wave 
only within the soliton, while it turns to the rest outside the wave, e.i. at 
C ^ — 00, and the only result of the action of the soliton is a shift of a particle 
by AC ~ A. 

The analysis of possible turbulent acceleration of dust particles by a random 
set of low-frequence Alfven solitons is out of the scape of our paper. One can 
speculate only that as far as translational heating of dust particles by low- 
frequency MHD waves is concerned, one can connect it with an ensemble of 
solitons passing randomly through a given volume of the interstellar space. 
In this picture the effective translational temperature is determined by the 
spectrum of low-frequency MHD solitons. 
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5 Summary 

In this paper we considered acceleration of charged dust particles by nonlinear 
low-frequency MHD waves under the conditions typical for interstellar plasma. 
We have shown that contrary to the case of intermediate frequencies Ud <C 
uj <^ uji described in [6], the acceleration efficiency is strongly limited by the 
amplitude of the nonlinear wave: while transversal velocities of dust particles 
Vdi_ are proportional to the amplitude of the soliton, the longitudinal velocities 
Vdz are of the second order of the amplitude. As a consequence, dust particles 
can be accelerated by such low-frequency waves to the velocities of few to tens 
percents of the Alfven speed. 
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conditions [1] and dissipation of magnetic flux [2] in the interstellar molecu- 
lar clouds. The dynamical importance of dust grains, and circulation of their 
mass (dust destruction and growth) in the interstellar medium is essentially 
determined by their kinetic temperature, and therefore understanding of how 
they are accelerated is of significant interest [31H] . In general, magnetohydro- 
dynamic (MHD) waves are known as an efficient heating source of plasma 
components. For instance, interaction of dispersive Alfven waves with plasma 
results in Joule heating of the electrons, generation of zonal flows and forma- 
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gyro- resonant acceleration [3] . More recently ponderomotive forces of the shear 
Alfven waves are shown to be a powerfull mechanism of dust acceleration [6], 
such that in a short time a microsized dust particle of mass ~ 10~^^ g can 
reach a few km s~^ velocity - comparable to the phase velocity of the Alfven 
wave. In [6] a steady nonuniform field envelope with a fixed field gradient 
(caused, for instance, by a nonuniform plasma density distribution along the 
wave vector) is described. Contrary, in our analyses below the ponderomotive 
force is unsteady in the rest reference frame where the pulse envelope moves 
with the group velocity. In this approach we are able to derive the perturbation 
amplitudes in the reference frame connected with the pulse, for densities and 
longitudinal velocities of plasma components driven by the mutual action of 
Coulomb and ponderomotive forces. Therefore, the nonuniformity of the field 
amplitude and the corresponding plasma acceleration in the pulse follows in 
our approach from a self-consistent solution of MHD equations. 

In [6] they considered a relatively high-frequency limit, corresponding to fre- 
quencies much higher than the dust cyclotron frequency Ud- It might be ex- 
pected though that in the low-frequency end u <ti ujd the presence of dust 
can have an important influence on dynamics of Alfven waves. In particu- 
lar, one can think that if dust particles are accelerated, the energy gained 
by them from the Alfven wave can be sufficient to decrease its amplitude, 
which in turn can reduce the acceleration efficiency. In this paper we explore 
the efficiency of acceleration of dust grains by the nonlinear low-frequency 
Alfven waves, and show that the nonlinearity strongly restricts the ability of 
Alfven waves to transfer their energy to charged dust grains. Such solitons 
can be actually treated as elementary objects of MHD turbulence in the low- 
frequency range. We follow therefore recent discussion of interaction of such 
low-frequency Alfven solitons with the plasma [7118] . 

In the conditions of the interstellar plasma this low-frequency range, uj -C 
~ 3 X 10~^^ Hz, can be of great interest from the point of view of dust ac- 
celeration because a considerable energy fraction of the interstellar turbulence 
can be locked in low- frequency motions [9]. Of course, nonlinear interactions 
move low-frequency Alfven perturbations toward higher frequences. However, 
it seems natural to assume that in the interstellar environment a steady state 
spectrum of perturbations establishes on a relatively short time scale (defi- 
nitely shorter than 100 Myr - the crossing time for a spiral wave), so that 
perturbations are always present in the low-frequency end. 

We proceed as follows. In Section 2 we find the dispersion relation and the 
group velocities for longitudinal Alfven and fast MHD waves, in Section 3 
we derive the nonlinear wave equation; Section 4 contains estimates of the 
velocities of dust grains accelerated by the Alfven and fast MHD soliton waves; 
summary of the results is given in Section 5. 
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In a fully ionized dusty plasma the dispersion relation for longitudinal Alfven 
and fast MHD modes can be written as 

iV' = ei ± 62, 



where N = kc/u, 
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a = e, z, d, LUpa and are the plasma and cyclotron frequencies of the ath 
component. In the low- frequency limit uj ujd equation ([T]) can be written in 
the form 
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where signs ± belong to the Alfven and fast MHD waves, respectively, VAi = 
Bo/^/4^Fpi is the Alfven speed. 

Assuming mdPdl ZdfriiPi ^ 1 and uj <^ pdOJdj P-, P = Pi + Pd, one can get from 
here 

= k\,l (l ± ^kVAd] , (4) 



where v\ = Bq/Atip, r^d = VA/^d- The group velocity and its derivative for 
Alfven waves are 

u = — = va \ l ^ TAdk , -j7 = —VArAd > 0, (5) 
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3 Nonlinear parabolic equation for low-frequency Alfven and fast 
MHD waves 



Following the procedure described in [TUfTT] one can derive the nonlinear 
Schrodinger equation for the envelope of the Alfven wave package in the form 

dEi u'd'^Ei u ^ ujlk^^ ( kv,k \ ^ 

+ 7717^2^^ — — hEi = 0, 7 



where 



U = ± ■ rADVA, 

Pi + Pd 



is the dispersion of the group velocity; sign "plus" belongs to the left-polarized, 
"minus" to the righ-polarized waves, Nk = ftk/nko is the normalized density 
perturbation of ki\i component. In the left-polarized wave ujk is positive for the 
electrons and dust particles, and negative for the ions, in the right-polarized 
Alfven mode is positive for the ions and negative for the electrons and dust 
particles. 

Assuming vti <^ t'A ^ one can neglect the electron inertia and write the 
equation of motion in the form 

dNe _ e d(j) ( d k d \ |Ep 



The equations of longitudinal motion of the ions and dust are 

dVz ^ e_d(^ ^ {—+ ^ ^\ |^|2 

dt M dz M'^uiui — uj) \dz uj uji — u dt j ' 
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Last terms in the r.h.s. of equations fl9l)- f|TTl) describe the ponderomotive force 
acting on the particles by the moving wave. 

With using expansion over a small parameter uj/ui -C 1 one can transform 
equation of motion of the ions and dust f|TOl) - f|TTl) in the reference frame moving 
with the wave group velocity to the form 

dV. edct^ f + F 
— u—— = — , 12 
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After integration of ^ and ffT^ - f[T^ one can find tlie density perturbations 

as 
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Then the quasineutrahty condition 

UioNi = UeoN^ + ZudoNd, (14) 



gives the following solution for the Coulomb potential 
e0 Wi W2 



and the explicit solutions for the density perturbations 



(15) 



Wo Wo ( TJM \ Wo 



here c^riioTe/neoM is the square of the ion-sound dust speed. The longitudinal 
velocities are connected through the continuity equations with the density 
perturbations by the equations 

v^k = uNk. (17) 



With accounting of equations (fT6l) and (fT7|l in the last terms in equation ([7]) 
one can reduce it to the following form for the new variable h = B/Bq = 
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The Lighthill criterion for the modulational instabihty fulfils when the signs of 
u' and of the nonlinear term in (fTSl) coincide. Thus, for Alfven wave with n' > 
modulational unstability takes place whenever the conditions c^pa/p < < 
hold. In the interstellar plasma with pd/p ~ 0.01 both these conditions are 
valid. For the right-polarized Alfven waves with u' < modulational instabil- 
ity can develop only in one of the two cases: either > cj. or < c^pd/p, 
which hardly can hold in the conditions of the interstellar matter. Therefore, 
left-polarized Alfven waves with the amplitude higher than a critical value 
experience in the interstellar plasma modulational instability which results in 
formation of a sequence of Alfven solitons. In the limit u ^ Cg the free term 
in (1181) becomes singular, which means that the quasineutrality condition is 
not valid anymore. 

Using the standard procedure we represent the magnetic field as & = fo_Le**, 
where b± = b±{( — vt), $ = $(C — ^i^), what means that the envelope solution 
we are seeking for drifts slowly (i.e. v <^ u) in the reference frame moving 
with the group velocity, while the phase $ drifts with the velocity vi ^ u. 
Separating the real and imaginary parts of (ITSll we arrive at 
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where 
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The solution of ffTOD is 



b± = b^ch-' (^^) , (21) 



where A = ^Ju'/ab'^, and vi can be found as vi = v{l — au'b'^/v'^)/2. In 
general, the solution (!2T]) is determined by the two free parametes: b^ and 
V. It is readily seen that collisions of dust particles with ions and atoms are 
unimportant on the scales of interest. Indeed, the characteristic length of the 



6 



soliton (I2T!) for uj <^ Ud can be estimated as ^ rAdb^, or for typical 
parameters in the interstellar plasma ^ lO^^b^ cm, which on the lower 
end is much less than the drag free path of dust particles: £ ~ 3 x lO^^n"^ 
cm, where the mass ratio rrid/m ~ 10^*^ and mean grain radius a = O.l/xm are 
assumed. 



4 Longitudinal and transversal acceleration of dust particles 



In the laboratory reference frame longitudinal acceleration of dust particles 
is determined by the last equation in (fT6|) and (fT7|) : \Vdz\ — u\b\'^/2, so that 
within the validity of the approximation the amplitude of dust velocity on 
nonlinear Alfven wave can be of several percents of the Alfven speed, which 
in the conditions of the interstellar plasma can vary from ~ 3 km s~^ in 
molecular gas (p ~ 10~^^ g cm~^, B ~ 10/uG) to ~ 10 km s~^ in diffuse 
HI phasei (p ~ 10~^^ g cm~^, B ~ 3yuG). Therefore, reasonable conservative 
estimate of dust acceleration along magnetic lines in such environments may 
be of 0.3 to 1 km s~^, respectively, which however depends on exact value of 
the soliton amplitude. This is much lower than the estimate obtained in [B]. 
Similarly, the dust velocity in perpendicular direction due to acceleration by 
the soliton can be estimated as \Vd±\ ^ u;\b±\/k u\b\, which is therefore a 
factor of larger than acceleration in the longitudinal direction. 

The acceleration of dust grains by a regular low-frequency Alfven soliton 
{uj ^ ujd) is connected with the ponderomotive force as in [6], and differs 
qualitatively from the mechanism described in [3], where translational heat- 
ing of dust particles is associated with stochastic acceleration of dust particles 
by MHD waves with ~ cj^ through cyclotron resonant interactions. Indeed, 
it is readily seen that dust particle moves under the action of soliton wave 
only within the soliton, while it turns to the rest outside the wave, e.i. at 
C ^ — 00, and the only result of the action of the soliton is a shift of a particle 
by AC ~ A. 

The analysis of possible turbulent acceleration of dust particles by a random 
set of low-frequence Alfven solitons is out of the scape of our paper. One can 
speculate only that as far as translational heating of dust particles by low- 
frequency MHD waves is concerned, one can connect it with an ensemble of 
solitons passing randomly through a given volume of the interstellar space. 
In this picture the effective translational temperature is determined by the 
spectrum of low-frequency MHD solitons. 
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5 Summary 

In this paper we considered acceleration of charged dust particles by nonlinear 
low-frequency MHD waves under the conditions typical for interstellar plasma. 
We have shown that contrary to the case of intermediate frequencies Ud <C 
uj <^ uji described in [6], the acceleration efficiency is strongly limited by the 
amplitude of the nonlinear wave: while transversal velocities of dust particles 
Vdi_ are proportional to the amplitude of the soliton, the longitudinal velocities 
Vdz are of the second order of the amplitude. As a consequence, dust particles 
can be accelerated by such low-frequency waves to the velocities of few to tens 
percents of the Alfven speed. 
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